Excised acoustic black holes: the scattering problem in the time domain 
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The scattering process of a dynamic perturbation impinging on a draining-tub model of an acoustic 
black hole is numerically solved in the time domain. Analogies with real black holes of General 
Relativity are explored by using recently developed mathematical tools involving finite elements 
methods, excision techniques, and constrained evolution schemes for strongly hyperbolic systems. In 
particular it is shown that superradiant scattering of a quasi-monochromatic wavepacket can produce 
strong amplification of the signal, offering the possibility of a significant extraction of rotational 
energy at suitable values of the angular frequency of the vortex and of the central frequency of the 
wavepacket. The results show that theoretical tools recently developed for gravitational waves can 
be brought to fruition in the study of other problems in which strong anisotropics are present. 
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I. INTRODUCTION 

General Relativity (GR) deals with dynamical de- 
formations of space and time, making a massive use 
of tensor calculus in parallel with continuum mechan- 
ics. In recent years the exploration of Einstein's the- 
ory has received a strong acceleration from the growing 
efforts in high-precision experiments involving gravity. 
Ground-based and space-born interferometers Hke LIGO, 
VIRGO, TAMA, GEO600, and LISA should at long last 
allow detection of gravitational waves. The analysis of 
the expected signals will require precise theoretical tem- 
plates to be compared with the measured waveforms via 
pattern matching. For this reason a variety of theo- 
retically simulated scenarios are needed, as obtained by 
numerical integration of the systems of partial differen- 
tial equations (PDE's) of GR. The striking developments 
of Numerical Relativity over the last ten years, which 
rely on presently available large-scale computational re- 
sources, have shown that some engineering tools such as 
the finite-elements techniques can be fruitfully adopted 
to these ends. On the other hand, the need to deal with 
event horizons, Cauchy horizons, and infinite-curvature 
singularities during the numerical evolution still poses 
serious theoretical and numerical challenges. 

In the past a very simple and yet profound analogy 
between GR and Newtonian physics has been noticed 
Q 1^ Q and it can be hoped that the former may benefit 
from the wide body of knowledge which is available for 
the latter. In particular, given a perfect barotropic and 
irrotational Newtonian fluid, perturbations of the veloc- 
ity potential with respect to a background solution have 
been analysed and shown to satisfy a second-order linear 
hyperbolic equation with non-constant coefficients. By 
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algebraic manipulations this wave equation can be rewrit- 
ten as a Klein-Gordon field propagating on a pseudo- 
Euclidean four-dimensional Riemannian manifold: an in- 
duced effective gravity is at work in the fluid, in which 
the speed of sound plays the role of the speed of light. 

GR, however, prescribes the existence of black holes. 
Surprisingly, some of the induced geometries in the ana- 
log formulation have the same structure as those of a 
black hole. Sound waves crossing the horizon of these 
acoustic counterparts will never come back out. At this 
level a practical problem arises. A black hole acts on 
light propagation as an anisotropic medium having an 
infinite refractive index at the horizon, and from studies 
of strongly anisotropic media, one can anticipate that a 
numerical study of this problem will face serious difhcul- 
ties. One may also ask whether black-hole behavior could 
emerge in other fields where hyperbolic wave equations 
are met, for instance in dealing with electromagnetism 
of anisotropic media or with perturbations of elastic ob- 
jects. 

Fortunately, Numerical Relativity has developed the 
proper instruments to deal with this type of problems, 
through the use of excision techniques together with 
a constrained evolution scheme for strongly hyperbolic 
problems. The main purpose of the present work is to 
use this methodology in the study of the simplest acous- 
tic black hole, i.e. the draining bathtub model of Visser 
While much can be borrowed from the extensive lit- 
erature in Numerical Relativity, the task of developing 
the apparatus needed to pursue this black-hole analog 
on fully quantitative grounds by numerical simulations 
does not seem to have been previously addressed. The 
relevance of these studies to superradiant extraction of 
rotational energy from a vortex in a Bose-Einstein con- 
densate has been discussed elsewhere 4] . 

The plan of the paper is as follows. In Sect. II we 
present a detailed analytical study of the draining bath- 
tub geometry in analogy with standard treatments of a 
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rotating Kerr black hole. In Sect. Ill we recast the hy- 
perbolic equation to a constrained symmetric system of 
first-order PDE's and solve a numerical scattering prob- 
lem on the basis of a finite-element method. In Sect. 
IV we analyse the time evolution of both compact and 
quasi- monochromatic sound-wave pulses, report on tests 
of the validity of the proposed procedure, and pay spe- 
cial attention to the scattering of a quasi-monochromatic 
wavepacket in the superradiant regime. Finally, Sect. V 
concludes the paper with a discussion of the main re- 
sults and an illustration of future perspectives regarding 
vortices in superfluids. 



II. THE MODEL 

We aim at investigating sound-wave scattering pro- 
cesses from a sonic black hole by studying the time evolu- 
tion of a scalar field in the presence of a vortex in the fluid 
flow. It is well known P, Q that sound propagation in 
such a fluid is described by a wave equation which, under 
the assumption of long- wavelength perturbations, can be 
mapped onto a Klein-Gordon equation associated with 
an effective relativistic curved space-time background (an 
acoustic metric). We focus our attention on the acoustic 
metric associated to the draining-bathtub model intro- 
duced by Visser |^ for a rotating acoustic black hole. 

The model is based on irrotational barotropic and 
incompressible Euler's equations describing a (2 -I- 1)- 
dimensional flow with a sink at the origin. The flow is 
taken as vorticity-free (apart from a possible ^-function 
contribution at the vortex core) and angular momentum 
is conserved. These constraints imply that the density 
of the fluid, the background pressure, and the speed of 
sound c are constant throughout the flow. The back- 
ground velocity potential must therefore have the form 



^/)(r, (j)) = — calog(r/a) + ^la^cf), 



(1) 



where a is a length scale associated with the "radius"" 
of the vortex horizon and 17 is the (constant) rotation 
frequency. 

Such a velocity potential, being discontinuous on going 
through 27r radians, is a multivalued function. Therefore 
it must be interpreted as being defined patch-wise on 
overlapping regions around the vortex core at r = 0. The 
velocity of the fluid is then given by 



v{r, (p) — V?/'(r, <p) = 



(2) 



The acoustic metric associated to this configuration is 



a2c2 + a4f^2 
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-dtdr (3) 



dr'' + 



It is easily shown that the metric (O exhibits an acoustic 
event horizon located at r = a, where the radial compo- 
nent of the fluid velocity exceeds the speed of sound, and 



an ergosphere at r = rerg = a-\/l + W'o? /c^ . A detailed 
analysis is given in Appendix IXI 

Linear perturbations i/'^^-' = 'f of the velocity potential 
satisfy the massless Klein-Gordon scalar wave equation 
VV^'I' = on this background, i.e. 
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This PDE is strongly hyperbolic, and its numerical in- 
tegration is best performed by resorting to the tools of 
Differential Geometry (and thus of General Relativity). 

Equation |0J is completely separable and we anal- 
yse it in the frequency domain by setting ^E* — 
^-iujt^im<t>^ikz p^^-^^ where k and the integer m are the ax- 
ial and azimuthal wavenumbers, respectively. The result 
is a second-order ordinary differential equation (ODE): 
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P = 0. 



Here m is taken as integer in order to avoid polydromy 
problems, while fc is a real number fixed by the boundary 
conditions along the z axis. It may be interesting to 
notice that by the Frobenius method for ODE's it can 
be shown that the indicial equation computed at r = 
a has two roots, Ai = and A2 = ^ {uj — mO). The 
latter expression leads to the well known condition for 
rotational superradiance: < a; < mfl for lo € TZ, as 
shown by the asymptotic study of the ODE for r = a 
and r ^ 00 HIEIl- 

Time dependent modes with m — can be analytically 
solved in terms of hypergeometric functions. In this spe- 
cial case Eq. © can be cast in a Schroedinger-like form, 
with a "dressing" potential outside the horizon. Following 
a standard procedure 01 , by rescaling the field as 



P(r) = 5(r)exp<i-i 
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and by introducing the "tortoise- like "coordinate 
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which maps the radial coordinate from [a, -f-oo] to 
[—00, +00], we obtain the equation 
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Figure 1: The horizon-dressing potential V{r)a^ is plotted as 
a function of r/a for ka — 0.02, 0.2, 0.4, and 0.6 (from bottom 
to top). 



where 

, ,2 5 1 /l + 4Pa2\ 3 
V{r) fc - T -jT - T ^ + o -T- 9) 

It should be noted that r, is a monotonic function of r 
and differs appreciably from r only in the vicinity of the 
horizon. 

In Fig. ^the potential V{r) is plotted as a function 
of r for several values of k. Waves of a given frequency 
u; are scattered by the real potential V{r) and can only 
partially tunnel towards the horizon for fc < 0.5. As fc 
increases further, the potential becomes more and more 
attractive and the back-reflecting barrier disappears. 

In the general case m ^ the equation is complex, and 
a simple Schroedinger-like picture no longer applies. Sev- 
eral numerical studies of the wave equation have re- 
cently been performed in the frequency domain for k = 
[^Q. However, these studies are based on a radial equa- 
tion which differs from Eq. ©. New time and angle 
coordinates have been introduced in order to minimize 
the dragging effects and to work in the region outside the 
horizon only. The physical interpretation of the results 
is not straightforward as the azimuthal quantum number 
associated with the new angular coordinate reflects dif- 
ferent geometrical properties from those of the cylindrical 
angle (j). This is a well-known problem in Kerr black hole 
physics lil . In order to avoid complications with the an- 
gular coordinate, we shall keep the original cylindrical 
coordinate system involving horizon-penetrating coordi- 
nates. A perturbative formalism using these coordinates 
is highly appropriate in connection with numerical rela- 
tivity simulations using black-hole excision. 



III. THEORETICAL FORMULATION OF THE 
SCATTERING PROCESS 

We solve Eq. 10} numerically in the time domain by 
resorting to some of the mathematical prescriptions that 



have been developed by Scheel et al. [lOj for the numer- 
ical integration of massless scalar-field perturbations on 
a rotating Kerr black hole background. 

The starting point is to introduce conjugate quantities 
of the scalar field 5*: 

where — (acr~^, — a^fJ, O) and a = c (see Appendix 
EI). These are used in Eq. to obtain a symmetric 
first-order hyperbolic system. We next factorize the de- 
pendence on z and 4> by setting 4" = ■0i(t, r)e""'^e'''^, 
n = 7^l(^,^)e™'^e^*=^ $i = 0i(t, r)e™'^e*'=^ $2 = im^ , 
and <I>3 = ik"^ . These transformations allow us to reduce 
the computational demands. The wave equation 10} is 
then replaced by a first-order set of coupled PDE's, 

dti^i + cdr (^1 — flTTi/r) — (ac — ima^ri) TTi/r^ -|- 

dfipi — cdr {0'4'i/^) ^ i'^c — ima^iVj ipi/^'^ ~ ctti (11) 
dt£,i + cdr (tti — aS^i/r) = 211710^^11^1 /r^ — ima^flS^i/r^ 

associated with the constraint 

|c| = |a,^i-(/.il = o (12) 

following from the definition of 4>i. Deviations of this 
quantity from zero provide a useful indicator of the qual- 
ity of the numerical results. 

Following the standard prescription for scattering pro- 
cesses in Kerr black holes Q, we choose as initial data a 
Gaussian pulse centered at r = tq and modulated by a 
monochromatic wave, 

i/)i(0,r) = Aexp[-(r ~ tq + ctf /b^ - ia(r - tq + ct)/c]\t=o 

(13) 

together with the associated initial conditions for tti 
and 4>i. The corresponding power spectrum is P{ijj) — 
-PmaxCxp [— (w — cr)262/4c2]. Thc modulated pulse thus 
has a Gaussian frequency distribution centered aX lo = a 
and the choice of a in the range < cr < mU, with m ^ 1 
corresponds to the superradiant regime. 

It is well known that the numerical integration of wave 
equations in cylindrical geometry must handle reflected 
waves close to the inner boundary (see e.g. Koyama tl l| ) . 
Similar difficulties arise in dealing with black holes, where 
theoretically correct boundary conditions are frequently 
in need of ad hoc modifications such as "background sub- 
traction" and "peeling off prop erties of the field (l/r" 
fall-offs)"at large distance [ij. The best way to solve 
the problem is to use penetrating coordinates or, better 
yet, the Kerr ingoing coordinates in Eq. I)A6|I . These 
permit to extend the field integration inside the horizon 
via the so-called excision technique. Inside the horizon, 
the light cones point towards the singularity, so that the 
interior region is causally disconnected from the region 
outside the horizon. Thanks to this property, one can 
impose a generic well-behaved boundary condition inside 
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the horizon and sidestep the problem of constraint vi- 
olations inside the acoustic black hole. This procedure 
is computationally very demanding, as it relies upon a 
very fine discretization around the sonic horizon. Ow- 
ing to requirements of high resolution, long integration 
times, and distant outer boundary, even the simple Klein- 
Gordon equation in Kerr black-hole geometry has only 
recently been successfully implemented with a powerful 
multi-processor code . Following Scheel et al. 0| , we 
have implemented the excision technique for an acoustic 
black hole and chosen the boundary conditions as follows. 

We have set no boundary conditions on the inner 
boundary r-m , with < rin < a. Regarding the outer 
boundary Tout, the proper criterion is to adopt condi- 
tions which will not violate constraints, will not sponta- 
neously generate unphysical waves, and will not signifi- 
cantly back-reflect any incoming signal. In our case we 
have first defined the directional derivatives of 'J along 
the principal null directions given in Eq. (|A5|I . namely 

u+ = n^a^^-cx (7rl-|-0l)e™'^e*'=^ 

u- = ;^9^^oc(^i-0i)e™'^e''=" . (14) 

At large distances we must impose a purely outgoing 
boundary condition, in order to avoid backscattered ra- 
diation which would contaminate the signal received at 
an observational point located closer to the vortex: con- 
sequently we set u~ ~ 0, which leads to no boundary 
conditions for tpi and to 0i = tti. In the limit of a flat 
space-time (a = 0), from Eq. (|10|l this condition becomes 
cdj-tpi +dt'ipi = and yields ipi = /i2(?' — ci), which is the 
familiar condition for zero ingoing (left-moving) waves. 
These boundary conditions did not lead to constraint vi- 
olations nor generation of spurious signals or of signifi- 
cant reflected waves during the numerical integration of 
the field equations. Moreover, by stopping the integra- 
tion before any physical signal has causally reached the 
outer boundary, even tiny reflected waves are completely 
eliminated. 



IV. NUMERICAL RESULTS 

The integration of the field equations is performed by 
using finite elements techniques instead of the standard 
finite-difference ones. We have used in particular the 
nonlinear engine of FEMLAB©, adopting a fifth-order 
polynomial Lagrange element with a uniform mesh in- 
terval with Ar = 0.05. Although this mesh is not very 
fine, the use of higher-order elements yields very good 
results. The integration is performed for r G [0, 150a] 
while time is in the range t € [0, 150a/c] with an op- 
timal time step chosen by the nonlinear solver. In or- 
der to avoid constraint violations and to secure stability, 
we have used finite elements of lower (third or even sec- 
ond) order close to the outer boundary. In our numerical 
analysis we have mainly focused on the scattering of a 
"quasi-monochromatic" (large b) wavepacket both inside 
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Figure 2: Constraint violations at the inner (a) and outer (b) 
boundary for m = in the case of a compact wavepacket 
(b = 0.5a, a — 0) for — lAa/c and ro — 10a. 

and outside the superradiant regime, although the case 
of a spatially "quasi-localised" (compact) wavepacket has 
also been examined. 

In the calculations we have adopted values of the pa- 
rameters that can be related to realistic superfluid sys- 
tems and examined both the non-superradiant and 
the superradiant regime by suitable choices of the value 
of a/fl. The angular frequency f2 was analysed in the 
range 0.14 < fia/c < 14, and the values m = and 1, 
ka = 0.02, and A — 0.3c were selected together with 
To = 10a and b — 0.5a for the quasi-localised packet and 
To = 50a and b = 10a for the quasi-monochromatic one. 
Code units corresponding to a = 1 and c = 1 have been 
used. 

As an example of our diagnostics. Fig Inland Fig El plot 
the evolution of |C| as a function of t for m ~ and 
TO = 1 at the inner and outer boundaries. These figures 
show that our boundary conditions preserve the constri- 
ants to a very satisfactory degree. In both cases the per- 
turbation does not grow indefinitely in time: the acous- 
tic black hole appears to be stable against these pertur- 
bations. The evolution of |C| for m = 1 at the outer 
boundary confirms that more attention has to be paid 
to possible constraint violations in this case (as already 
found in However, such violations in Fig. |3Ib) are 

still negligible. In both cases, constraint violations re- 
main within the range 10~^ to 10"^", which means that 
causality violations, if present, are practically irrelevant. 

As a further illustration of the stability of our results, 
we compare in Fig. ^the values for the time-dependent 
energy of the wavepacket calculated with the standard 
mesh of 0.05 with those obtained on the finer mesh with 
Ar = 0.003. There evidently is complete agreement be- 
tween the two sets of results, even in the time range of the 
scattering process where the energy of the wavepacket is 
varying very rapidly in time. 

We can now describe the main physical results of the 
numerical integration. Fig[5land Fig|Sl report snapshots 
of the time evolution of the real part of ipi as a function of 
r/a for tc/a = 0, 5, 7 and 10 in the case to = and m — 1, 
respectively. The peak in the signal at t = shows the 
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Figure 3: Constraint violations at the inner (a) and outer (b) 
boundary for m — 1. The parameters used in the calculation 
are as in Fig. |21 



ingoing Gaussian pulse. At t = 5a /c the pulse starts 
being affected by the dressing potential and the radiation 
is backscattered near t — 7 ioi m = and t = 10 for 
m = 1 . The small reflected wave then relaxes towards the 
steady state (Re^'i = 0) after the characteristic ringing- 
mode oscillations. 

Figure [71 shows a density plot of the real part of ipi in 
the range r/a ^ [0, 60] and tc/a £ [0, 60] for the case m = 
(no superradiance, left) and m = 1 (superradiance, 
right). 

The initial Gaussian pulse moves towards the vortex 
horizon placed at r = a, its trajectory being bent by 
the potential outside the horizon. The bending of the 
pulse trajectory, with light cones heading towards the 
horizon, is consistent with similar findings in numerical 
relativity hj^ . Qualitatively similar results are obtained 
in the non-superradiant regime for m = 1 by choosing 
cr = 1.2n. 

Finally, the superradiant behavior of the acoustic black 
hole during the scattering of the wavepacket is illus- 
trated in Fig. |S1 and El In FigjSl we show a typical 
time evolution of the energy of wavepacket Ep{t) — 



rH 



(pM/2)/p" d0/(f dz/J^^%^dr with VI = V0(i), p 
the density and H the axial extent of the fluid, nor- 
malized to its initial value Ep{0), for a = 0.7c/a and 
= \Ac/a, within the superradiant regime (m = 1, 
solid line) and outside it (m = 0, dashed fine). In 
the non-superradiant case, the energy of the scattered 
wavepacket goes asymptotically to zero, indicating that 
all the impinging energy is lost to the vortex sink. In 
the superradiant case instead, the energy of the back- 
scattered wavepacket exceeds its initial value, showing 
extraction of energy from the ergosphere at the expense 
of the rotational kinetic energy. The energy gained via 
superradiance is seen to exceed in this case twenty per- 
cent of the initial value Ep(0). 

In Fig. [51 we show the superradiance efficiency, de- 
fined as the ratio between the total energy gain AEp — 
(i?p(oo) — Ep{0)) of the wavepacket and the background 
energy E^,, for three values of a/fl in the superradiant 
range, as a function of fla/c. In the perturbative regime 
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Figure 4: The energy Ep{t) as a function of tc/a, normalised 
to its initial value Ep{0), as calculated with two different 
choices of the mesh (see text). The parameters used in the cal- 
culation are a = 0.07c/a, Q — O.lAc/a, b — 10a, and ro = 50a. 
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Figure 5: Snapshots of the time evolution of the perturbation 
'i/'i representing a compact wavepacket at t = (solid line), 
t — 5 (dotted line), t = 7 (dashed line), and t = 10 (dot- 
dashed line), as a function of the distance r from the the 
vortex center for the case m = (t in units of a/c). The solid 
line shows the initial Gaussian pulse travelling towards the 
vortex. The parameters used in the calculation are as in Fig. 



(for U < 3c/ a, say) the efficiency of energy extraction 
grows very rapidly with fl, especially at large values of 
the ratio cr/ri between the central frequency of the quasi- 
monochromatic wavepacket and the angular rotation fre- 
quency of the vortex. Although substantial values of the 
efficiency are beyond the scope of the perturbative Klein- 
Gordon description used in this work, the indication from 
Fig. 1^1 is that an efficiency approaching unity may be 
obtained by a suitable match between the frequency of 
the wavepacket and the rotational frequency. Indeed, it 
seems that nonlinear effects may primarily determine the 
way in which the energy extraction behaves as it becomes 
comparable to the background energy. 
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Figure 6; Snapshots of the time evolution of the perturbation 
-i/ji representing a compact wavepacket at t = (solid line), 
t = 5 (dotted line), t = 7 (dashed line), and t — 10 (dot- 
dashed line), as a function of the distance r from the the 
vortex center for the case m = 1 (t in units of a/c). The 
parameters used in the calculation are as in Fig. |21 
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Figure 7: A density plot of the real part of i/'i in the r-t 
plane for m — (left) and m — 1 (right) in the case of a 
quasi-monochromatic wavepacket. The calculations refer to 
n = 3c/ a and a = 0.711 (color online). 



V. DISCUSSION 

Interest in pursuing analogue models of gravitational 
physics has been steadily growing in the recent past. Our 
attention in the present work has been directed to sound 
wave propagation in a background fluid flow as an ana- 
log for field propagation in a curved space-time, as first 
highlighted in the seminal work of Unruh . As a first 
step in a long-term program, we have studied the scat- 
tering and radiance phenomena from a black hole whose 
background space-time may be modelled by the flow-field 
configuration associated with a vortex excitation in a 



fluid with a drain. We have deliberately chosen values 
of the system parameters that are typical of superfluids 
in a Bose-Einstein condensed state. The inverse crossing 
time of a vortex in such a dilute gaseous superfluid is 
cja K, 15KHz, of the same order as for a cosmic black 
hole of radius a w 10km. Such a close match is a con- 
sequence of the very low speed of sound in condensates, 
which is a few mm/s. 

Our analysis of the acoustic analog has shown that 
acoustic black holes behave under several respects like 
black holes in GR. They have an event horizon and an 
ergosphere, and can superradiate at the expense of their 
rotational kinetic energy. Once excited, they come back 
to their initial state and are stable against sharp pertur- 
bations. The special methods that one typically uses in 
the theoretical treatment of black holes in space-time can 
be applied here in a straightforward way. Our calcula- 
tions have indicated that substantial superradiance effi- 
ciencies may be achieved from a "terrestrial black hole" by 
a proper choice of the system parameters. This sugges- 
tion will have to be substantiated by the inclusion of 
nonlinear and quantum effects, and indeed the acoustic 
analogy is still awaiting for a deeper analysis and under- 
standing. 

An especially important direction for further applica- 
tions of numerical relativity to superfluids will be to de- 
velop a proper inclusion of vortex quantization. As is 
well known, the vorticity of a superfluid is quantized in 
units of the circulation k — h/Al, i.e. fin = nh/{Ma?), 
where n is an integer, M the particle mass, and a the 
size of the vortex core. Consistency requires that an 
integer number of quanta be exchanged between a vor- 
tqx and a wavepacket in a scattering process. A semi- 
classical treatment is thus justified only in the case of 
scattering from so-called giant vortices built from up to 
n ~ 60 quanta, as are indeed observed in some experi- 
ments on Bose-Einstein condensates inside strongly an- 
harmonic confinements (see e.g. It will also be 

necessary in this context to consider a vortex with no 
drain, such a configuration being accessible in the lin- 
ear regime by an immediate modification of the present 
acoustic analog 0,^^. More generally, the inclusion of 
both nonlinear and quantum effects in wavepacket prop- 
agation inside superfluid flows of dilute quantum gases 
can be based on the use of the time-dependent Gross- 
Pitaevskii equation. Advanced techniques are available 
for its numerical solution in a variety of configurations 
[TtI I. Work along these lines is in progress. 
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Figure 8: Time evolution of the energy gain Ep{t) of the 
wavepacket (in units of its initial value _Ep(0)) for a — 0.7c/a 
and SI = lAc/a, for the superradiant case (m = 1, solid line) 
and the non-superradiant case (m = 0, dashed line). The 
time t is in units of a/c. 




Figure 9: Efficiency of superradiance, defined as the ratio 
between AEp = _Ep(oo) — Ep{0) and the background energy 
Eh, plotted on a semi-logarithmic scale as a function of the 
angular frequency Q (in units of c/a) for three values of the 
ratio a/Q. 



Appendix A: PROPERTIES OF THE ACOUSTIC 
METRIC IN EQUATION (3) 

The line element ^ is most conveniently analysed by 
recasting it in the following (3 + l)-form: 

ds^ = -a^dt^ + 7y {dx' + P'dt) {dx^ + (3^dt) . (Al) 

The three-metric is the flat Euclidean space in cylindrical 
coordinates 7^ = Diag(l, t-'^, 1). Latin letters are used 
for spatial indices that are raised and lowered with jij 
and its inverse 7*-' = Diag(l, r~^, 1), whereas Greek let- 
ters indicate four-dimensional quantities. The shift is 
given by — {acr^^ , —a^ftr"^ ,0), and consequently 
Pi = {acr~^, —a^n,0). The lapse is a = c. 

The slices t = const are flat, and consequently this 
effective space-time is written in Painleve— GuUstrand- 
like form. The presence of a shift vector manifests the 



well know phenomenon of frame dragging. The metric 
depends on the radius only, so that three Killing vectors 
are associated to cyclic coordinates. This, together with 
Group Theory, ensures that both geodesies (sound rays) 
and wave fields can be solved by separation of variables. 
Moreover this Killing vector is not rotation free, i.e. the 
space-time is stationary. The metric, which diverges in 
the origin only, possesses an event horizon, and in order 
to demonstrate this assertion we have to analyse the null 
geodesies. 

We start by writing the principal function 
S{t,r,cj),z) = -Et + f{r) + L(j) + Qz, where the 
quantities E,L, and Q are constants of the motion. The 
associated Hamilton- Jacobi equation iVSVpS* = for 
massless particles becomes 



1 - 



dr 



2a {Er^ - a^VLL) df 



E^r^ - 2a^nLEr^ 



dr 



- (A2) 

= 



and is solvable by quadrature in terms of elliptic inte- 
grals. The gradient of S leads to the quantity [/^ — V^S*, 
which is conventionally called the four- velocity of the null 
geodesies, although for light-rays it should more appro- 
priately be called the four dimensional wave vector k^. 

If we assume no motion along the z axis, i.e. Q = 
0, and zero angular momentum L = (at infinity the 
motion is purely radial) , Eq. ljA3p above delivers the two 
solutions 



/i,2(r) = T—r 
c 



Ea 
c 



The associated four velocities are 



rE 



E Eila^ 



(r + a) ' c ' rc^ (r -|- a) 
rE E EQa^ 



,0 



(A3) 



, (A4) 



(r — a) ' c ' rc^ [r — a) 



,0 



The quantity n'^ is the vector tangent to an ingoing con- 
gruence of null geodesies, and the quantity l^^ corresponds 
instead to the outgoing ones. The ingoing null signals are 
regular on the surface r = a, whereas the outgoing ones 
diverge. This means that at r = a every signal can en- 
ter but not escape, i.e. there is an event horizon. This 
pair of congruences of null vectors are the principal null 
directions associated with this geometry. More in detail, 
they automatically bring the Weyl tensor into its canon- 
ical form, i.e. a Petrov type-D manifold as for station- 
ary black holes in General Relativity. However, unlike 
GR cases, the Ricci scalar is non-zero, implying that the 
Goldberg-Sachs theorem for type-D space-times does not 
apply. 

A new coordinate system for this acoustic geometry 
can be tailored to the ingoing light rays associated with 
these vector fields. This step is easily performed by defin- 
ing a new set of coordinates v and ■0, which are related 
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to the old ones by 



V r a , ( r\ 

+ - log 1 + - , 

c c c V a) 

iP log — — 

c V ^ + a 



(A5) 



and the ingoing congruence becomes n'' = (0, —E,0,0), 
as expected. The space-time is now expressed in "Kcrr- 
ingoing like" coordinates. 



Geometrically, the new coordinates compress time and 
untwist the angle in the neighborhood of the horizon. 
The corresponding metric simplifies to 



a2 a4^2 
1 - ^ 



dw^ + 2dt;dr (A6) 



2na' 



dvdip + r^dV^^ + dz^ 



The vanishing of the gu component of the metric tensor 
means that the norm of the temporal Killing vector = 



{■§f)'^ changes sign at the radius re 
This manifests the presence of an ergosphere and sug- 
gests the possibility of rotational energy extraction via 
superradiance. 
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